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Recent elastic and inelastic neutron scattering studies of the highly frustrated pyrochlore antifer- 
romagnet Tb2Ti207 have shown some very intriguing features that cannot be modeled by the local 
(111) classical Ising model, naively expected to describe this system at low temperatures. Using the 
random phase approximation to take into account fluctuations between the ground state doublet 
and the first excited doublet, we successfully describe the elastic neutron scattering pattern and 
dispersion relations in Tb2Ti207, semi-quantitatively consistent with experimental observations. 
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The search for the physical realization of a spin liq- 
uid in three dimensions has been a long quest for the 
condensed matter community. Recently, geometrically 
frustrated magnetic systems have been the focus of in- 
tensive experimental and theoretical studies because it 
is believed that geometrical frustration can inhibit the 
formation of long-range order, thus enabling the sys- 
tem to remain paramagnetic down to low temperatures. 
Among three-dimensional systems, the pyrochlore lat- 
tice of corner-sharing tetrahedra have been studied ex- 
tensively. It has been shown theoretically and numer- 
ically 0, 0, that for classical Heisenberg spins with 
nearest-neighbor antiferromagnetic (AF) exchange, there 
is no transition to long-range magnetic order at finite 
temperature. This makes AF materials based on the 
pyrochlore lattice excellent candidates to search for the 
low-tcmpcrature spin liquid state. A number of experi- 
mental studies on insulating pyrochlore materials have 
been carried out. Interestingly, most materials either 
develop long-range Neel order, such as Gd2Ti 2 07 Q 
and E^Ti^Cv or reveal spin-glass behavior, such as 
Y 2 Mo 2 7 m The "spin ice" materials, Ho 2 Ti 2 7 0,0 
and Dy 2 Ti 2 07 0, Il0| exhibit low-temperature thermo- 
dynamic properties reminiscent of Pauling's "water ice 
model" 11]. In these systems, an effective ferromagnetic 
interaction is frustrated due to the single-ion (111) Ising 
anisotropy 0, 0] . The behavior of these systems can 
be quantitatively described by the (111) Ising spin model 
with nearest-neighbor exchange and long-range dipolar 
interactions HElEH- 

Tb 2 Ti 2 07 shows, however, very different and intrigu- 
ing behavior. It is believed that Tb 2 Ti 2 C>7 belongs to 
the same family of (111) Ising systems as Dy 2 Ti 2 07 and 
Ho 2 Ti 2 07 but with an effective nearest neighbor AF in- 
teraction HHHIllIUm. The same s P in modc l that 
very successfully described the spin-ice systems predicts 
it to have a noncollinear Neel Q = order, with all 
spins pointing into or out of each tetrahedron, at about 
1 K |l0j . In dramatic contrast, Tb 2 Ti 2 07 remains a 
spin liquid, or "cooperative paramagnet" , down to 70mK 



|13l llfij. In addition, recent paramagnetic neutron scat- 
tering studies show that the scattering pattern for this 
material is not consistent with a (111) Ising model, while 
a Heisenberg spin model 0, ^} or some level of relax- 
ation away from the (111) Ising model 18] can better 
describe the observed neutron scattering pattern. These 
results suggest that the restoration of spin isotropy in 
the system, despite its expected Ising-like nature at low 
temperature 0, 0, 1^3, El > is essential in understand- 
ing the paramagnetic spin correlations. Inelastic neutron 
scattering studies have also been performed on this sys- 
tem and partial softening of the magnetic excitations at 
an energy of about 20 K has been observed [151 Hfl Eij ] . 
This has been attributed to a (spin) roton-like mode, as 
in liquid 4 He 0, 0, |23| , which further indicates a more 
isotropic nature of the spins. Given the ensemble of evi- 
dences, it would appear that one needs a more isotropic 
spin model to understand the paramagnetic spin correla- 
tions in this system. More importantly, such a "restora- 
tion" of spin isotropy may also be the key to understand- 
ing why Tb 2 Ti 2 07 fails to order down to low tempera- 
tures. In this paper , w e employ the random-phase ap- 
proximation (RPA) [24j to take into account the fluc- 
tuations between the ground state doublet and the first 
excited doublet. We successfully describe the observed 
paramagnetic spin correlations in Tb 2 Ti 2 07, without any 
assumptions regarding the nature of the spins, while we 
still obtain the Q = Neel order at low temperatures. 
This result makes the fact that Tb 2 Ti 2 07 fails to order 
even more puzzling. 

We begin with the model spin Hamiltonian, 

n = \ J2 s*Vaf (m^ + E^fM), (1) 

where i,j are indices of the Bravais lattice vectors for the 
FCC lattice, a, b are indices of sublattice basis vectors 
and a,/3 are indices for the spatial coordinates. TCcf is 
the single-ion crystal field (CF) Hamiltonian. The spin- 
spin interaction matrix J , including both exchange and 
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long-range dipolar interactions, reads 
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where 8 nn refers to nearest-neighbor interaction, and 
R"j' a denotes the a component of the interspin vector 
RJj that connects spins Si j£l and Sj.b. The nearest- 
neighbor exchange for Tb2Ti207 is given by J\ ~ 
—0.167 K [25| . and D^d = (.9Ms) 2 Mo/4vr, where /xg is 
the Bohr magneton, /iq the magnetic permeability, and 
5 = 3/2 for Tb 3 +. 

The single-ion susceptibility is given by [3| 
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where n v is the thermal occupation fraction for state v. 
The matrix elements for the single-ion states are given 
by M "^,a = EsWS^X.a, wh ere u% a is the rotation 
matrix from the local (a) frame defined on sublattice a to 
the global (a) frame. The spin operator S a acts on the 
CF states defined along the local quantization axis (l4| . 
The local wave function structure has been obtained from 
both experimental measurements and theoretical calcula- 
tions [1J, |2(J, |2jJ . The primary components are J z eigen- 
states I ±4) for the ground state doublet and | ±5} for the 
first-excited state doublet. The doublets are separated by 
an anisotropy gap A m 20 K, which is comparable with 
the Curie- Weiss temperature #cw ~ — 19 K [l4j . 

The RPA equation which takes into account the two- 
ion interaction contribution to the full susceptibility [24| 
is given by 

+ E x° Q ' Q7 HJa 7 /(q)x^(q^) = <Ux°< a », 

7,<5,c 

where J{<£) is the Fourier transformation of the interac- 
tion matrix J(i,f). The slowly converging, infinite lat- 
tice sum of the dipolar interaction is handled using Ewald 
summation techniques 0)13 • We solve for x"f(q, ui) 
numerically using LAPACK routines. 

The elastic and dynamical neutron cross-sections are 
related to the spin susceptibilities by summing over the 
sublattice contributions and taking the transverse com- 
ponents of S perpendicular to Q 24], 
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FIG. 1: (a) Experimental elastic neutron scattering pat- 
tern of Tb2Ti2C>7 in the (h,h,l) plane of reciprocal space at 
T=9 K, from Ref. Hfil . Dark blue shows the lowest intensity 
level, red-brown the highest, (b) S(Q) for the (111) Ising 
model at T=9 K. (c) S(Q) for the doublet-doublet model at 
T=9 K. Note that magnetic form factor |F(Q)| 2 is divided 
out in both experimental and theoretical results. 
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where Q = G + q, G is a reciprocal lattice vector of the 
FCC lattice, and q is a primitive vector in the first Bril- 
louin zone. r a is the sublattice basis vector and F(Q) is 
the magnetic form factor for Tb 3+ . 

Figure d shows the experimental scattering pattern 
in the (h,h,l) plane at 9 K 0|. Note that there is a 
strong intensity maximum around (0,0,2). Figured 3 
shows the S(Q) calculation using the (111) Ising model, 
i.e., the anisotropy gap A — > oo. It is clear that this 
model fails to reproduce the correct neutron scatter- 
ing pattern as observed experimentally, as there is only 
weak non-critical intensity around (0,0,2). It has been 
shown that in the (111) Ising model, the intensity at 
(0, 0, 2) and (0, 0, 0) are exactly correlated, and vanishes 
for Tb 2 Ti 2 7 0, EH E3 • Using a more realistic doublet- 
doublet model for Tb2Ti207, we are able to qualitatively 
reproduce the experimentally observed scattering pattern 
(Fig- EO ■ It captures most details of the experimental 
pattern, such as the intensity maximum around (0, 0, 2), 
the minima around (0,0,0), (2,2,2) and (0,0,4). This 
excellent agreement between theory and experiment in- 
dicates that to properly understand the spin correlations 
in Tb 2 Ti 2 07, the fluctuations out of the ground state 
doublet to the first excited doublet, originating from the 
first term in Eq. J5J), are important. It should be noted 
that this result is obtained from the simple level scheme 
described above with no other assumptions about the na- 
ture of the spin or the details of the wave functions. 

Line scans in reciprocal space along three high- 
symmetry directions (0,0, 1), (h,h,0) and (h,h,h) at 
T =4 K are plotted in Fig. [2J Open symbols are the 
data points for S(Q). Filled symbols are the dispersion 
relations E{Q) for the lowest-energy band of magnetic 
excitations at 4 K. The dispersion relations show good 
semi-quantitative agreement with the experimental dis- 
persion observed in single crystal Tb 2 Ti 2 0y reported in 
Fig. 9 of Ref.0. We also find that the minimum in E(Q) 
corresponds to the maximum in 5(Q) in all symmetry 
directions, as is observed in the experiment. The sharp 
jumps in dispersions, e.g., near (0,0, 1) and (0,0,3), cor- 
respond to shifts of spectral maximum between different 
branches, since the experiment and current calculation 
track only the maximum peaks in the energy spectra. 
In principle, there should be four branches of magnetic 
excitations with different intensities due to the four sub- 
lattice structure. Higher resolution experiments need to 
be performed to map out the dispersion relations of these 
branches. 

To study the temperature effects on the dispersion re- 
lations of the lowest three bands of magnetic excitations, 
we use the wave functions obtained from crystal field cal- 
culations in the point charge approximation and in the 
fixed J = 6 manifold for Tb 3+ 14]. This gives the CF 
level scheme with two low-lying doublets as described 
above, and two higher-energy singlets The en- 

ergy levels are given as 0, 24.6, 113.9 and 132. 7K, con- 
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FIG. 2: The open symbols show the line scans of the S(Q) 
for the three high-symmetry directions: from top to bottom, 
(0, 0, 1), (h,h,0) and (h,h,h). The filled symbols show the 
dispersion relation for the lowest-lying branch of magnetic 
excitations at T = 4 K. 



140 r 



12K 



30K 



135- 



130 ( 
120 



115 

110, 

30 

25 
20, 






1 


2 


3 


4 


1 1 1 1 1 1 1 
1,1,1, 





1 


2 


3 


4 






^^—^ 


w 







i 


2 


3 


4 



(hhh) 



FIG. 3: The dispersion relations of the three lowest-lying 
magnetic excitations along the (h, h, h) direction at both 12 K 
and 30 K. The lowest band displays a more pronounced en- 
ergy dispersion as temperature is lowered, in agreement with 
experimental observations [l6j|. 



sistent with experiments. The two higher singlets have 
large | ± 3) contributions and some mixing with other 
| J, Mj) components [l4|. All the excited states are con- 
nected to the ground state doublet through the S + and 
S~ operations and the excitations are visible via neutron 
spectroscopy • Figure shows the dispersion rela- 
tions of the three lowest-lying magnetic excitations along 
(h, h, h) at 12 K and 30 K at the same energy scale. The 
lowest band becomes more dispersive as the temperature 
is lowered, while the higher two bands do not show much 
change with temperature. 
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Fluctuations between the doublets are important in 
understanding the experimentally observed scattering 
pattern. These fluctuations originate in the first term 
of Eq. , where the anisotropy gap A = Ei — E enters 
the formulation algebraically instead of exponentially as 
in the elastic (second) term through the thermal occupa- 
tion fraction. If the two-ion interaction, characterized by 
Oqw, is comparable with A, then the effect of this term 
is large and part of the isotropic response is restored. By 
setting A — > oo or the matrix element between two dou- 
blets to zero, this term is eliminated, and the scatte ring 
pattern is reduced to that of the (111) Ising model [FT!. 
This explains why the (111) Ising model works so well 
in describing the spin correlations in spin ice materials 
HE3, since in these systems A « 200 ~ 300 K [U H3 
and 6»cw ~ 0.5 - 2 K (jEHOU, so the effects of the fluc- 
tuations out of the ground state doublet are negligible. 
To understand the different temperature dependencies of 
the dispersion for the three lowest-lying magnetic exci- 
tations, we inspect the thermal factor n„ M = n„ — n M in 
the first term of Eq. (J2J). For temperatures much lower 
than A, only the ground state doublet is populated, and 
this factor is almost the same for all three bands. When 
the temperature reaches a value T « A, the first excited 
state doublet begins to be populated, therefore reduc- 
ing the factor uqi = tiq — m, while those for the two 
higher levels remain essentially unchanged. This change 
in thermal occupation factor results in a large tempera- 
ture dependence of the dispersion in the lowest band. 

The low temperature ground state of our model in RPA 
is the noncollinear "all-in" or "all-out" Q = Neel state 
[l(i| with T c w 1.8 K. Exact T c value depends on the 
details of the crystal field wave functions. Suprisingly, 
this value is very close to the T c obtained recently in 
hydrostatic pressure measurements for pressure greater 
than 1 GPa [23|. By standard Holstein-Primakoff ex- 
pansion aroud this ground state we obtain the spin- wave 
excitation spectrum with a large gap and the reduction 
of the staggered magnetization due to quantum fluctu- 
ations Am/m < 10 _4 [2^|. This result should be taken 
only as an indicator, though, since the ground state dou- 
blet (Mj = ±4) is not at its saturation value Mj = ±6 
due to the crystal field interaction. The magnetic ground 
state doublet is a non-Kramers doublet, and the states 
are time-conjugate of each other, so there exist no matrix 
elements between them [2flj | , indicating that any quantum 
fluctuations causing spins to flip result from a higher- 
order virtual process via excited crystal field levels. 

In conclusion, with a simple level scheme and the RPA 
to include the fluctuations between the ground state 
and the first excited state doublets, we are able to de- 
scribe semi-quantitatively the experimentally observed 
neutron scattering pattern and the energy dispersion in 
Tb2Ti2C>7. Our results indicate that the crystal field ef- 
fects are important in this system due to the fact that the 
anisotropy gap and the two-ion interaction are compara- 



ble, and that some isotropy is present in the response. 
The ground state of our model, however, is still the "all- 
in" or "all-out" Q = Neel state with T c w 1.8 K, 
which is stable against the lowest order quantum fluctua- 
tions. These results indicate that a mechanism that may 
restore more isotropy and thus increase the quantum fluc- 
tuations leading to the suppression of T c , requires further 
detailed study. 
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